I 


AD-A263  999 


NPS-MA-93-012 

NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


NONLINEAR  MODE  COUPLING 
IN  FREE  ELECTRON  LASERS 
by 

C.  L.  Frenzen 

Technical  Report  For  Period 
January  1993  -  March  1993 

Approved  for  public  release;  distribution  unlimited 

Prepared  for:  Naval  Postgraduate  School 
Monterey,  CA  93943 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


Rear  Admiral  T.A.  Mercer  Harrison  Shull 

Superintendent  Provost 

This  report  was  prepared  in  conjunction  with  research  conducted  for  the  Naval  Postgraduate 
School  and  funded  by  the  Naval  Postgraduate  School. 

Reproduction  of  all  or  part  of  this  report  is  authorized. 

This  report  was  prepared  by: 

r.L 

C.L.  FRENZEN 
Associate  Professor 


Reviewed  by:  Released  by: 


*  unclassified 


SECURIT-  CLASSiF'CATiON  Of  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


form  Approvrtf 
0MB  r^o  0704  0I8S 


1«  REPORT  SECURITY  ClASSlPlCATlON 

unclassified 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER{S) 


lb  RESTRICTIVE  MARK  NGS 


3  Distribution/ ava.labil. TV  of  REPORT 

Approved  for  public  release, 
Distribution  unlimited 


5  MONITORING  organization  REPORT  NuMBERiSi 


NPS-MA-93-012 


6a  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School 


6c.  ADDRESS  (Ofy,  Stare,  and  ZIP  Code) 

Monterey,  CA  939A3 


NPS-MA-93-012 


6b  OFFICE  SYMBOL  7a  NAME  OF  MONiTORiNG  ORGAN  ZATiON 
(If  app/icabfe) 

MA  Naval  Postgraduate  School 


7b  ADDRESS  (Ofy,  State,  and /IP  Code; 

Monterey,  CA  939A3 


8a  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Naval  Postgraduate  School 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Monterey,  CA  939A3 


8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicibit) 


O.M  &  N 


'0  SOURCE  OF  funD'NG  Numbers 


PROGRAM 
element  no 


ACCES5iO^  NO 


1 1  TITLE  (Include  Security  Classification) 

Nonlinear  Mode  Coupling  in  Free  Electron  Lasers 


12  PERSONAL  AUTHOR(S) 

C.L.  Frenzen 


13a  type  of  REPORT 

Technical 


16  supplementary  notation 


13b  TIME  COVERED 

'=«OM_i-93 _ TO  3-a3 


14  DATE  Of  REPORT  (Year.  Month  Day) 

93  March  11  _ 


CCSATI  CODES 


GROUP  SUB-GROUP 


18  SUBJECT  TERMS  (Confmue  on  reverie  tf  neceii^ry  identity  by  b^ock  number) 

Free  electron  lasers,  optical  wave 


19  abstract  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  optical  field  in  free  electron  lasers  can  sometimes  develop  sidebands  and  exhibit 
verr’  complex  behavior.  The  purpose  of  this  report  is  to  develop  a  two  mode  model  of  the 
optical  field  in  a  free  electron  laser  in  whidi  the  optical  modes  are  coupled  by  the 
free  electrons  in  the  laser.  This  model  may  be  used  to  study  the  long  term  stability 
of  the  optical  wave  and  the  developnent  of  sidebands  in  the  optical  field.  In  this 
report  the  equations  of  motion  governing  the  interaction  of  an  electron  with  the  optical 
modes  will  be  derived. 


20  DISTRIBUTION/ availability  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

09  UNCLASSIFIEO/UNLIMITED  □  SAME  AS  RPT  □  DTIC  USERS  I  UnclSSSified 


OD  Form  1473,  JUN  86 


Previous  editions  are  obsolete 

S/N  0102-LF-0U-6603 


SECuRlTv  CLASSIFICATION  Of  This  =aGE 


Nonlinear  Mode  Coupling  in  Free  Electron 

Lasers 

C.L.  Frenzen 

Department  of  Mathematics 
Naval  Postgraduate  School 
Monterey  CA,  93943 

Abstract 

The  optical  field  in  free  electron  lasers  can  sometimes  develop  side¬ 
bands  and  exhibit  very  complex  behavior.  The  purpose  of  this  report 
is  to  develop  a  two  mode  model  of  the  optical  field  in  a  free  electron 
laser  in  which  the  optical  modes  are  coupled  by  the  free  electrons  in 
the  laser.  This  model  may  be  used  to  study  the  long  term  stability 
of  the  optical  wave  and  the  development  of  sidebands  in  the  optical 
field.  In  this  report  the  equations  of  motion  governing  the  interaction 
of  an  electron  with  the  optical  modes  will  be  derived. 

1.  Introduction 

Initial  theoretical  studies  of  the  free  electron  laser  (FEL)  used  quantum  me¬ 
chanics  ([l]-[3]).  However  it  was  showm  in  [4]  that  a  classical  treatment  of 
FELs  would  adequately  describe  most  phases  of  operation  ([5]-[10])  and  re¬ 
sult  in  substantial  simplification  of  the  resulting  equations  describing  their 
operation. 

One  way  to  understand  some  of  the  main  features  of  the  FEL  problem 
is  to  consider  the  motion  of  a  pendulum  with  weak  damping  and  a  nearly 
constant  tangential  forcing.  In  a  mathematical  model  of  the  FEL  with  a 
tapered  magnetic  wiggler,  however,  the  weak  damping  term  in  the  pendulum 
equation  cannot  be  ignored  because  the  wiggler  is  extremely  long  compared  to 
its  period,  and  one  must  account  for  the  “cumulative”  effects  of  this  damping 


over  the  entire  length  of  the  wiggler  in  order  to  obtain  accurate  solutions. 
The  only  way  to  systematically  account  for  these  higher  order  effects  is  to 
formulate  a  completely  dimensionless  problem  which  contains  dimensionless 
paurameters  whose  magnitudes  characterize  various  distinguished  limits. 

In  a  FEL  frequently  the  optical  field  can  develop  sidebands  and  become 
highly  complex  in  its  behavior.  The  purpose  of  this  report  is  to  develop  a 
model  in  which  there  are  only  two  modes  present  in  the  optical  field:  these 
modes  are  coupled  to  one  another  through  the  electrons  in  the  FEL.  This 
model  can  be  used  to  study  the  long  term  stability  of  the  optical  wave  and  the 
development  of  sidebands  in  the  optical  field.  A  regime  of  particular  interest 
is  that  of  sustained  resonance  in  the  two  mode  model.  In  this  report  we  shall 
derive  the  equations  of  motion  governing  the  interaction  of  an  electron  with 
the  optical  fields  in  the  two  mode  model  of  a  FEL. 

2.1  Equations  of  Motion 

We  shall  formulate  the  equations  of  motion  as  a  first  order  system  of 
ordinary  differential  equations.  Our  derivation  follows  that  in  [11].  However 
we  shall  be  extending  the  derivation  given  there  for  a  single  optical  mode 
to  the  two  optical  mode  case.  The  convention  we  shall  use  is  that  upper¬ 
case  letters  will  represent  dimensional  variables  and  lower-case  letters  will 
represent  dimensionless  variables. 

To  determine  the  leading  order  equations,  consider  a  single  electron  mov¬ 
ing  in  a  helical  wiggler  magnetic  field  given  by 

B,„  =  5u,{xcos^Z-l-Ysin|^z},  (1) 

where  By,  is  the  magnetic  field  amplitude  and  is  its  wavelength.  The 
vectors  X,  Y,and  Z  are  orthonormal  and  are  associated  with  the  coordinates 
X,  Y,  and  Z,  where  Z  is  the  coordinate  along  the  direction  of  the  magnet  axis. 
We  shall  assume  that  the  radiation  (or  signal)  field  is  a  linear  superposition 
of  two  modes  and  is  given  by 

E.  =  Ei‘>  +  Ei'\  (2) 

where 

E<;>  =  -Ei'HZ)  {Xsin  e.  -I-  Y  cos  0,}  ,  »•  =  1, 2  (3) 
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and 

B.  =  B<’>  +  Bi^), 

(4) 

where 

Bi‘>  =  Z  X  =  1,2. 

(5) 

is  the  radiation  amplitude  of  mode  i  and 

0.  =  |^z-u»<‘)r  +  ^<;),i  =  i,2 

(6) 

where 

is  the  wavelength,  is  the  frequency,  T  is  time,  and 

is  the 

phase  of  mode  i,  i  =  1,2. 

W2  shall  assume  that  the  electron  is  relativistic,  though  its  transverse 
velocity  (the  velocity  orthogonal  to  the  Z  aj<is)  is  small.  Let  (3  be  the  dimen¬ 
sionless  electron  velocity  (/3  is  the  electron  velocity  divided  by  c,  the  speed 
of  light).  We  shall  decompose  (3  into  a  component  parallel  and  a  component 
perpendicular  to  the  Z  axis: 

(3  =^||Z  +  ^x,  (7) 

The  equation  of  motion  (Newton’s  second  law)  for  the  electron  is 
P  =  eE,-\-el3  x(Ba-fBu,)+space  charge  force-f  radiation  reaction  force,  (8) 


where  e  is  the  electron  charge  and  P  =  jmc^  is  the  electron  momentum.  7 
is  the  Lorentz  factor,  defined  by 

c  is  the  speed  of  light,  and  m  is  the  electron  rest  mass.  Criteria  for  neglecting 
the  effects  of  the  radiation  reaction  force  and  the  space  charge  force  were  given 
in  [12],  and  were  derived  by  examining  higher  order  terms  in  the  equations 
of  motion.  These  conditions  will  be  given  later  in  this  report. 

The  conservation  of  energy  implies  (7mc’)  =  ec^  •  E,  or 

(10) 

To  gain  a  rough  idea  of  the  approach  to  be  followed,  we  shall  assume 


DD 


shall  repeat  the  same  calculations  in  more  detail  with  two  modes  in  the 
radiation  field. 

The  momentum  equation  (8)  implies 

(7mc/3  )  =  7{mc/3  )  +  7mc/3  =  7^(mc;3  )(/3  •  /3  )  +  7mc/3  ,  (11) 

since  7  =  7^y3  •  4  . 

Combining  the  momentum  equation  (11)  and  the  energy  equation  (10), 
we  conclude  that 

7mc4  =  eE,  —  e/3  (/9  •  E,)  +  e/3  x  (Z  x  E,)  +  e/3  x  (12) 

Since  ^  X  (Z  X  E,)  =  Z(/3  •  F,)  —  E,(/3  •  Z),  equation  (12)  may  be  written 
as 

7mc/3  =  e(l  — /3||)E, -t- cZ/3  •  E,  +  e/3  x  Bu,.  (13) 

Using  the  fact  that  1  —  ;5||  1  and  that 

/3  ■E.=/3x-E.  (14) 

we  obtain  the  approximate  equation 

/3  =  — /3  X  Bu,.  (15) 

7mc 


Substituting  (3  ,  the  solution  of  equation  (15),  into  the  energy  equation  (10), 
we  find 


eE,/3x  .  ^ 

7  =  - sin  <t>, 

(16) 

me 

•  27rc  A„,  1  +  6^^ 

A„  i‘-A,  2y  /■ 

(17) 

where  <f>  —  (2jr/ +  27r/A,)Z  —  and 

“  2irmc2  ’ 

(18) 

is  defined  as  the  dimensionless  amplitude  of  the  wiggler  in  Section  (2.2).  The 
details  of  the  above  calculations  (for  the  case  of  two  modes)  will  be  given  in 
Section  (2.3). 
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We  can  now  define  a  (dimensionless)  resonant  electron  energy  7r  for  which 
0  =  0.  This  energy  is  obtained  from  equation  (17)  and  has  the  value 


T'r  = 


■^u;  1  + 


A. 


(19) 


A  free  electron  laser  usually  operates  in  a  state  for  which  7  as  7,  (see  [2],  [6], 
and  [7]).  Equation  (16)  indicates  that  electrons  with  (dimensionless)  energy 
7r  in  a  magnetic  wiggler  with  wavelength  and  field  amplitude  will 
spontaneously  emit  radiation  (if  >  <  0)  or  absorb  radiation  (if  >  >  0)  of 
wavelength 


A,= 


A.(l-}-62) 


27? 


(20) 


2.2  Normalization  and  Nondimensionalization 


It  is  very  important  to  decide  on  characteristic  scales  for  the  dependent 
and  independent  variables  in  the  study  of  a  system  of  differential  equations 
by  asymptotic  methods.  The  appropriate  choice  of  scales  is  usually  dictated 
by  the  physical  regime  of  interest  represented  in  the  system  under  study  and 
can,  in  fact,  be  changed  to  suit  a  different  regime.  They  key  point  is  that 
once  a  set  of  scales  is  chosen  all  dependent  and  independent  variables  become 
dimensionless  and  order  one  in  magnitude.  There  then  is  no  ambiguity  about 
the  size  or  numerical  magnitude  of  the  dependent  and  independent  variables. 
Of  course,  there  are  now  dimensionless  parameters  present  in  the  problem 
which  are  functions  of  the  characteristic  scales.  It  is  f/ieir  orders  of  magnitude 
which  indicate  what  physical  regime  is  under  study,  and  what  is  small,  etc. 

We  shall  use  X^l27r  =  To  as  a  characteristic  unit  of  length  and  To/c  =  Tq 
as  a  characteristic  unit  of  time.  Let  2irfXu,  =  /^u,  be  the  wavenumber  of  the 
wiggler.  Then  in  terms  of  these  units  we  define  the  dimensionless  variables 


z  —  K^Z  = 


t  =  K^cT  = 


To 


(21) 


(Note  that  lower  case  variables  will  always  be  dimensionless.) 
We  now  introduce  the  dimensionless  field  variables: 


K^mc^ 


B,,  Cj  — 


K^mc^ 


(22) 
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where  mc^  =  0.511  MeV.  Using  the  normalizations  introduced  in  equations 
(21)  and  (22),  equations  (8)  and  (10)  may  be  written  in  dimensionless  form 
as 

.)  =  e,  +  X  (b,  +  b„,)  +  higher  order  terms,  (23) 

at 

^  =  (24) 

In  equation  (23),  the  higher  order  terms  include  the  normalized  space  charge 
force  and  the  normalized  radiation  force. 

To  characterize  the  order  of  magnitude  of  each  term  in  these  equations 
we  need  to  introduce  a  small  parameter  £,0  <  e  1.  One  way  to  do  this 
is  to  choose  the  parameter  £  to  be  a  measure  of  the  slow  variation  of  the 
wiggler  field  along  the  magnet  axis.  To  do  this  we  define  £  as 


£  = 


2nL 


U 

— ,  or  £ 


1 

/ 


(25) 


where  L  and  1  are  respectively  the  dimensional  and  dimensionless  length  of  the 
wiggler  and  Lq  is  the  characteristic  length  defined  by  the  wiggler  in  equation 
(21).  The  smallness  of  £  signifies  that  the  magnetic  wiggler  produces  a  small 
effect  which,  for  an  electron  moving  along  the  magnet  axis,  lasts  for  a  long 
time. 

Using  the  small  parameter  £  we  can  define  a  “slow”  variable  z  as 


z  =  ez  ,  0<£'Cl. 


(26) 


The  slow  variable  z  enables  us  to  define  “slowly  varying”  functions  of  z  in  the 
following  sense:  an  0(1)  function  of  z  which  varies  slowly  along  the  wiggler 
axis  while  maintaining  its  order  of  magnitude  can  be  written  as  f(r).  Note 
that  df/dz  =  0{e).  A  slowly  varying  function  f  of  the  dimensionless  variable 
z  may  be  written  ^ls 

/(i)  =  /(£^)  =  /(^£)  =  /(|), 

so  that  slowly  varying  means  changes  occur  over  length  scales  which  are  much 
smaller  than  the  length  of  the  wiggler. 

Typically  the  wiggler  length  L  is  several  meters  and  the  wavelength  A,^ 
of  the  wiggler’s  magnetic  field  is  on  the  order  of  a  centimeter,  so  £  is  a  small 
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number,  about  10“^.  Also  since  the  dimensionless  wiggler  length  1  is  !/£,  the 
wiggler  must  be  considered  long  in  the  mathematical  sense.  This  means  that 
the  cumulative  effects  of  higher-order  terms  may  be  important. 

We  will  now  denote  all  normalized  (dimensionless)  variables  with  an  over¬ 
bar.  The  wiggler  magnetic  field  is  then  written  in  normalized  form  as 

b„  =  5,„(2){xcosz -!-y  sinr}.  (27) 

It  can  be  seen  now  from  equation  (19)  that  b^u  I  leads  to  7^  ie., 

extra  high  energy  electrons  are  required.  On  the  other  hand,  if  6,^  1, 

the  wiggler  field  is  then  too  weak  to  significantly  decelerate  electrons.  We 
conclude  that  the  only  reasonable  choice  is  =  ^>^(2)  =  0(1).  Because  the 
operation  of  free  electron  lasers  is  near  the  resonant  state  defined  in  equation 
(19)  (see  references  [2],  [6]  ,  [7]),  and  since  typical  values  of  7r  are  large, 
we  shall  write  the  dimensionless  energy  of  an  electron  and  the  dimensionless 
wavenumber  of  the  signal  in  the  following  form  consistent  with  equation  (19): 

7  =  (28) 

Y  =  (29) 

In  the  abtpve  equation,  fc,  is  the  0(1)  dimensionless  signal  wavenumber 
and  £“  =  /z  is  a  second  small  parameter  which  measures  the  ratio  X,fX^. 
The  determination  of  fi  comes  from  a  consideration  of  the  physical  regime 
under  discussion.  Recall  that  variables  with  an  overbar  are  dimensionless 
and  0(1). 

After  this  preliminary  analysis,  let  us  assume  again  that  there  are  two 
modes  present  in  the  radiation  field.  In  normalized  form  its  components 
appear  as 


6^9  =  — =  Ixsin^j +ycos^i), 

(30) 

b<'>  =  z  X  ef), 

(31) 

1,  2,  and  where 

ft  =  e-“Si'>(z  -  ()  + 

(32) 

again  for  i  =  1,2.  Note  that  the  total  radiation  field  is  given,  as  before,  by 
equations  (2)  and  (4).  The  parameter  (  in  equation  (30)  is  again  determined 
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by  the  physical  regime  under  discussion.  We  shall  show  in  section  [2.3]  that 
the  only  really  interesting  limit  for  the  free  electron  laser  is  ^  =  0. 

The  radiation  parameters  and  depend  on  the  wave  evolution  pro¬ 
cess  and  must  be  determined  by  Maxwell’s  equations.  There  are  some  open 
questions  here  which  will  be  discussed  later.  We  will  see  that  e,  evolves  on 
a  slower  scale  than  b^. 

There  are  basically  three  types  of  magnetic  wigglers  or  undulators.  There 
are  those  in  which  varies,  those  in  which  varies,  and  those  in  which 
both  and  b^  vary.  Our  approach  to  this  problem  can  be  adapted  to  any 
of  these  wigglers.  However,  the  work  in  [7]  suggests  that  the  constant  Au, 
case  is  probably  the  best  method  for  enhancing  efficiency.  Since  it  also  offers 
apparent  hardware  advantages  we  will  restrict  our  attention  to  this  case. 

2.3  Higher  Order  Correction  Terms 

In  this  section  we  calculate  higher  order  correction  terms  for  the  equations 
of  motion  and  we  mention  conditions  derived  in  [12]  for  neglecting  the  space- 
charge  force  and  the  radiation  reaction  force. 

We  begin  by  looking  at  equation  governing  the  temporal  evolution  of  the 
perpendicular  component  of  the  velocity,  From  equations  (23)  and  (24) 
we  find 

-I-  (^it^  +  0±)  X  (i  X  e#  bw).  (33) 

Since  eg  is  perpendicular  to  z,  we  find  that  the  last  equation  may  be 
written 


e  “^^^(/3x)  =  (1  -  0\\)^s  -t-  i^|Z  X  bw  +  (  higher  order  terms)  (34) 

Now  ,  for  a  single  electron,  ^  =  )9j|.  We  can  regard  this  as  an  ordinary 
differential  equation  expressing  z  as  a  function  of  t.  A  little  later  it  will  be 
convenient  to  invert  this  relationship  and  to  express  t  as  a  function  of  z. 

From  [12],  the  following  conditions  must  hold  in  order  to  neglect  the 
radiation  force: 

=  (35) 

Conditions  necesssary  for  neglecting  the  space  charge  force  are 
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(36) 


Once  Q  is  determined  (and  it  will  be  shown  later  that  (  is  zero),  the  above 
two  equations  provide  conditions  which  must  hold  in  order  to  consistently 
neglect  the  effects  of  the  radiation  and  space  charge  forces  on  an  electron.  In 
the  second  equation,  the  ratio  ^  is  the  ratio  of  the  plasma  frequency  to  the 
signal  frequency.  It  will  turn  out  that  this  ratio  is  a  fundamental  parameter 
measuring  the  rate  of  evolution  of  the  radiation  field.  Consequently  the 
second  condition  on  the  neglect  of  the  radiation  force  pla3’s  m  important 
role  in  the  evolution  of  the  optical  field. 

Inserting  the  specific  field  quantities  into  equation  (34)  then  gives 

c”"''^^(/3x)  =  -/?j|6,^(z)(5-sinz  -  ycosz) 

-(1  -  (x  sin^i  +  ycos^i)  (37) 

sin  ^2  +  ycos$2) 

+higher  order  terms, 
where  6^  =  since  bi* ^  =  i  x  el* ^ 

Noting  that  dt  =  dz,  we  integrate  equation  (37)  to  obtain  the  form 

(f  cos  2  +  y  sin  z) 

7 

t' 

— t^^^^(xsin2  —  ycosz) 

7 

— ^xcos^i  —  y  sin^i)  (38) 

“fki 

1(2) 

_^3o/2+€_^  ^2  -  y  sin  ^2) 

7k» 

where  6,  =  e“"fc6*(2  —  t)  +  In  equation  (38),  b'^  stands  for  the  derivative 
of  b^  with  respect  to  €°z.  However  since  is  a  function  of  2  =  cz,  it  follows 


> 
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that 


(39) 


d(e°z) 

The  error  in  equation  (38)  is  due  to  ignoring  the  variations  of 

d6ti,/<f(e“*),6„  and  with  2.  To  obtain  (38)  ,  we  have  also  assumed  that 
W0-)  =  ie.,  the  transverse  velocity  of  the  incoming  beam  can  be 

ignored  before  it  enters  the  magnetic  wiggler. 

Now 

J  •  e,  =  /?x  •  (ei'^  +  e<2))  (40) 

since  the  component  of  the  electron  velocity  parallel  to  the  wiggler  axis  is 
orthogonal  to  the  signal  fields.  Inserting  equations  (37)  and  (38)  into  the 
above  expression  in  (40)  ,  we  find  after  some  calculation  and  rearrangement 
that 


%  I 

gO+t  ■  sin(2  -f  ^i)  +  sin(g  +  ^2) 


_  2a+{ 


bj:^ 


cos(r  +  ^1)  +  . cos(r  +  O2) 


rrl^2) 

^2a+2<  + 


sin(f>j  —  62) 


+0(e=^). 


We  shall  now  define  the  phases  i  =  1,2  by 

=  2  +  0”  =  z  +  -  <)  +  (42) 

We  observe  that  since  many  of  the  parameters  in  the  problem  depend  on 
z,  as  previously  mentioned,  it  will  be  convenient  to  use  d/dz  rather  than  d/d* 
for  the  equations  of  motion.  The  velocity  dz/dt  can  be  determined  from  the 
relationships  dz/dt  =  and 
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Using  the  above  equation  and  the  expressions  for  7  and  in  equations 
(9)  and  (38)  ,  we  find  after  some  effort,  that 


^11 


_c2‘>+« 


cos  4>1 


^  ‘‘-V  A 


+0(c^) 


(44) 


Since 


it  follows  that 


^7  dz  dy 
di  dt  dz 

(45) 

d^  1  d^ 
dz  ;3||  dt 

(46) 

We  thus  can  combine  equations  (41)  and  (44)  above  to  derive  an  expression 
for  the  rate  of  change  of  7  with  respect  to  the  dimensionless  distance  along 
the  wiggler  axis,  z.  We  find  after  some  calculation  that 


d'f 

dz 


-e  ^  1  — rr~  sm  <Pi  4 - sin  <f>2 


— c 


2a+< 


l+t 

27' 


feu;?  . 

- sin  (^1  + 

sin  4>2 

7 

7 

— C 


2o+< 


~  cos  4>i  +  cos  <^2^ 


+t 


2o+2« 


sin(t^2-<^.)f^-^ 


+(?(c'^+^). 


(47) 


Equation  (47)  determines  the  evolution  of  7  with  respect  to  z,  and  it  can 
be  seen  that  this  depends  on  the  field  quantities,  and  on  the  phases 
defined  in  (42).  We  now  determine  the  evolution  of  these  latter  quantities 
with  respect  to  z.  Differentiating  the  expression  for  <f>i  in  equation  (42)  with 
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respect  to  z,  recalling  that  ^  =  l/i3||,  and  using  equation  (44)  for  we 
find  that 


d4>i 

dz 


+0{£'“) 


We  can  obtain  an  expression  for  the  rate  of  change  of  with  respect  to  z  in 
a  similar  fashion: 


d<t>2 

dz 


yi_  cos  <?i,  +  cos  <;i2  +  (49) 

+0(c2^), 


where  —  d(f>i*^ldz  must  be  determined  by  using  Maxwell’s  equations  to 
determine  the  evolution  of  the  optical  field. 

Note  that  the  electron  energy  (7  in  (47)  depends  on  both  the  phases 
4>i>4>2,  whose  evolution  is  determined  by  equations  (48)  and  (49).  Note  also 
that  the  evolution  of  these  phases  depends  on  the  electron  energy  7. 

If  we  now  compare  the  higher  order  terms  in  the  exprv  "ions  for  d^fdz, 
d<^i/dz,  and  d(f>2/dz,  we  see  that  unless  ^  =  0,  oscillatory  solutions  for  7,<^i, 
and  d>2  will  not  exist.  Since  the  existence  of  oscillations  in  7  implies  a  transfer 
of  energy  from  the  electron  to  the  optical  field,  and  this  is  essential  to  the 
operation  of  the  FEL,  we  find  that  we  must  take  ^  =  0  in  equations  (47), 
(48),  and  (49). 

At  this  point  we  have  two  independent  small  parameters:  t  and  c°  =  p 
(as  we  have  not  yet  specified  a).  From  equation  (25)  c  is  a  meaisure  of  the 
smallness  of  Ai„  compared  to  L,  while  measures  the  energy  7  of  the 

electron  (see  (28)).  As  in  [12],  near  the  resonant  state  (7  w  7^),  there  are 
two  distinguished  limits  corresponding  to  q  =  0  and  a  =  2/3.  The  latter 
case  is  known  as  “slow  resonance”  and  is  substantially  more  difficult  to  solve 
than  the  case  when  o  =  1,  and  for  this  reason  we  shall  restrict  our  attention 
to  the  a  =  1  case. 
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Equations  (47),  (48),  and  (49)  form  three  coupled  ordinary  differntial 
equations  for  7,  (pi,  and  <i>2.  The  coefficients  in  these  DDEs  are  slowly  varying 
functions  of  5  =  ez.  Though  the  single  mode  case  (see  [12])  can  be  cast 
into  Hamiltonian  form,  it  does  not  seem  possible  to  express  the  two  mode 
equations  in  Hamiltonian  form.  One  must  also  study  how  the  optical  wave 
interacts  with  and  evolves  according  to  the  presence  of  the  electron  beam. 
(We  have  derived  the  equations  for  a  single  electron,  but  generally  in  the 
electron  beam  of  a  FEL  there  are  many  electrons). 

In  order  to  maintain  the  coherence  and  monochromicity  of  the  radiation, 
the  evolution  of  the  optical  wave  should  be  slow.  Thus  the  leading  order  parts 
of  the  signal  parameters  6^1  and  4>\  must  either  be  constant  or  depend  on 
the  prescribed  slowly  varying  wiggler  parameters.  These  signal  parameters 
in  the  radiation  field  are  determined  by  the  wave  evolution  process  which  is 
governed  by  Maxwell’s  equations.  Though  determination  of  the  asymptotic 
behavior  of  the  signal  parameters  is  relatively  straight  forward  in  the  one 
mode  model  (see  [12]  for  example),  it  still  is  not  quite  clear  how  the  wave 
evolution  occurs  in  the  two  mode  model. 

We  leave  the  further  details  of  this  analysis,  and  the  study  of  the  three 
ODEs  in  (47),  (48),  and  (49)  to  a  future  report. 
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